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ABSTRACT 

We present the first measurement of the evolution of the galaxy group stellar mass function (GrSMF) to 
redshift z > 1 and low masses (M* > 10 12 M Q ). Our results are based on early data from the Carnegie-Spitzer- 
IMACS (CSI) Survey, utilizing low-resolution spectra and broadband optical/near-IR photometry to measure 
redshifts for a 3.6/im selected sample of 37,000 galaxies over a 5.3 deg 2 area to z ~ 1 .2. Employing a standard 
friends-of-friends algorithm for all galaxies more massive than logM*/M Q = 10.5, we find a total of ~ 4000 
groups. Correcting for spectroscopic incompleteness (including slit collisions), we build cumulative stellar 
mass functions for these groups in redshift bins at z > 0.35, comparing to the z = and z > mass functions 
from various group and cluster samples. Our derived mass functions match up well with z > 0.35 X-ray 
selected clusters, and strong evolution is evident at all masses over the past 8 Gyr. Given the already low level 
of star formation activity in galaxies at these masses, we therefore attribute most of the observed growth in 
the GrSMF to group-group and group-galaxy mergers, in accordance with qualitative notions of hierarchical 
structure formation. Given the factor 3-10 increase in the number density of groups and clusters with M + > 
10 12 M© since z = 1 and the strong anticorrelation between star formation activity and environmental density, 
this late-time growth in group-sized halos may therefore be an important contributor to the structural and star- 
formation evolution of massive galaxies over the past 8 Gyr. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: groups: general — galaxies: 
high-redshift 



1. INTRODUCTION 

Just as the relative roles of genetics and environment in 
shaping human populations are hotly debated, so are the anal- 
ogous processes argued among those who study the evolu- 
tion of galaxies. The "nature" of a galaxy is thought to be 
primarily determined by the mass of its dark matter halo, 
since in isolation the collapse of a given amount of matter 
should give rise to descendents wi th a consistent set of prop- 
erties ( Berlind & Weinberg 2002). However, like humans, 
galaxies don't often mature in isolation: their "nurturing" is 
brought about through gas inflows from the extended cosmic 
web, mergers with other galaxies, and even direct interactions 
with hot, dense intergalactic gas (e.g. ram press ure stripping; 
iKoopmann & Kennevl I2004t ICrowl et all 12005b. While it is 
well-established that galaxies ' morphologies (IDresslerlll980t 
Wilm anetaLll2009t) . colors dCooper et all 120071: iPatel et al I 
2009^7 and star-formation properties dKauffmann et al.l l2004; 
Pate letani20Til lOuadri et al.ll20l2T) differ in overdense and 
underdense environments, the reasons behind these correla- 
tions are still uncertain. 

By targeting mass-selected galaxy samples in a range of 
environments, the role of "neighborhood" can be disentan- 
gled from that of mass. On the extremes, void galaxies which 
evolve essentially in isolation tend to exhibit stable, g as-rich, 
star-forming configurations (Grogin & Geller 2000), while 
those in rich clusters are overwhelmingly quiescent and have 
morphologies which suggest frequent high-speed interactions 
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or mergers (e.g., Far ouki & Shapiroll981l) . Although such en- 
vironments provide invaluable laboratories for the study of 
ex- and in-situ galaxy evolution respectively, they represent 
rare and extreme neighborhoods. Most of the galaxy popula- 
tion in fact resides in the proverbial suburbs: relatively small 
galaxy groups containing a nywhere betwee n two to a few 
dozen significant members (lEke et al.ll2.004t >. Groups' influ- 
ence on their member galaxies may be dramatic, like quench- 
ing sta r formation by keeping col d gas from reaching galaxies 
(e.g., Ivan den Bosch et all 120081) and/or triggering bursts of 
star forma tion through slower gal axy-galaxy interactions than 
clusters (Larson & Tinslevl ll978l) . In turn, the assembly of 
groups (and the quenching processes therein) since z = 1 could 
be a major factor in the declin e of the cosmic star formation 
rate o ver the same time period dLillv et al.l[l9 96; Mada uet al.l 

cm. 

Despite the ubiquity of groups, homogeneous samples are 
notoriously difficult to assemble beyond the nearby universe 
(z > 0.4). Since each group contains only a few galaxies 
spread over ~ 1 Mpc or more, standard cluster selection tech- 
niques (e.g. angular overdensities exhibiting a well-defined 
red sequence) cannot be applied. Comprehensive spectro- 
scopic surveys over wide areas are the only reliable way 
to weed out groups from interlopers; however, at high red- 
shifts optically-selected spectroscopic surveys are biased to- 
ward UV-bright, highly star-forming galaxies at the expense 
of the massive red galaxies which dominate dense environ- 
ments. An infrared-selected survey, effectively producing a 
mass-limited sample at z ~ 1, is therefore necessary to pro- 
duce an unbiased census of massive galaxies at these redshifts 
dvan Dokkum et al.ll2"0()6t iKelson etal .1120 121) . 

The Carnegie-Spitzer-IMACS (CSI) Surve>Q currently un- 

3 More information at |http : / /csi ■ obs ■ carnegiescience . edu| 
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derway at the Magellan-Baade 6.5m telescope in Chile, has 
been specifically designed to characterize massive galaxies 
and their environments up to z ~ 1 .4. Here we present a mea- 
surement of the evolution of the group stellar mass function 
(GrSMF), to lower group masses and higher redshifts than has 
been previously achieved, using early data from CSI. Cosmo- 
logical parameters h = 0.7, fi m = 0.3, and = 0.7 are as- 
sumed throughout. 

2. OBSERVATIONS AND SPECTRAL FITTING 

In a companion paper (Kelson et all l2012h . we describe 
the CSI observing setup and strategy, data reduction, and 
SED fitting; a brief summary follows. We employ a sim- 
ple 3.6^m flux limit of m\s < 21 to select galaxies from 
three of the Spitzer Wide-Area Extragalactic Survey (SWIRE) 
fields: ELAIS-S1, Chandra Deep Field South, and XMM- 
Large Scale Structure (XMM-LSS), excluding stars with sim- 
ple optical/IR color cuts. The 3.6fim band lies in the rest- 
frame near-IR for galaxies at z = 0.5 - 1 .5, resulting in an ap- 
proximately mass-selected sample unbiased by dust and star 
formation. The first two years of observations focused primar- 
ily on the XMM-LSS field, giving it the best spectroscopic 
and photometric coverage at present; here we focus only on 
data from this 5.3 deg 2 field in our analysis. 

Targets were observed in multi-slit, nod-and-shuffle mode 
dGlazebrook & Bland-Hawthornl I200ll) with a multiplexing 
efficiency of ~ 1800 objects per 28' diameter mask, using the 
low-dispersion prism (LDP) on the I namori-Magellan Ar eal 
Camera and Spectrograph (IMACS: Ibressler et ail 1201 ll) in 
2009. In 2010, a second pass was mostly completed with 
the new uniform-dispersion prism (UDP) which provides su- 
perior resolution beyond A > 7500A. Some overlap between 
the LDP- and UDP-observed samples was included for cross- 
calibration and testing; overall, two passes have provided 
~ 40% completeness. The IMACS spectra are supplemented 
with moderate-depth J and K imaging from the NOAO Ex- 
tremely Wide-Field InfraRed Mosaic (NEWFIRM) and opti- 
cal photometry from the public CFHT Legacy Survey. 

The spectra were optimally extracted and combined, and 
prism and broadband IR data were jointl y fit using a general - 
ized set of starburst models described bv lKelson et al.l (120121) . 
With a series of b enchmarks — spectr oscopic redshifts in 
this field (VVDS; iLeFevreet ail 120051) . galaxies observed 
with both the LDP and UDP, and compa ring prism redshifts 
of ph ysically associated galaxy pairs (cf. lOuadri & Williams 
2010) — we infer typical redshift uncertainties of Az/(1 + 
z) ~ 0.005-0.015 at z < 1. Notably, unlike photometric red- 
shifts, these errors are comparable for red and blue galax- 
ies since the prism spectrum pinpoints both the 4000A break 
and emission lines. After excluding objects with poor SED 
fits (typically due to strong AGN components, bad photome- 
try, and/or incorrect slit placement) the final catalog contains 
about 37,000 galaxy redshifts over 5.3 deg 2 , with an effective 
stellar mass limit of about M t >3x 10 10 M Q at z ~ 1. 

3. GALAXY GROUPS AT 0.5 < Z < 1.2 

3.1. Finding groups with friends-of -friends 

Groups typically appear as associations of two to a few 
tens of galaxies over an area of a few hundred projected 
kpc; accurate redshifts and large areas are therefore criti- 
cal to building robust gr oup samples. We emp loy a simple 
friends-of-friends (FoF; iHuchra & Geller|[l982l) algorithm to 
find groups in the CSI catalog, limiting the sample to galaxies 



with Mi, > 3.2 x 10 10 M Q , the CSI mass limit for red galax- 
ies at z ~ 1.2. Even at the low end our redshift uncertainty 
(~ 0.5%) is larger than the expected velocity dispersion of 
most virialized groups and clusters. Redshift uncertainties 
were estimated (using the techniques mentioned in Section|2]i 
as a function of r-band magnitude and the quadrature sum 
of each galaxy pair's ler uncertainties was adopted as the 
redshift linking length. We used a transverse link ing length 
of 1.0/(l + z)Mpc (matching the I Yang et ail 120071 length of 
1 Mpc for SDSS groups). The mean group redshift was also 
adjusted as new members were added, and previous group 
members falling outside 1 a- from the new mean were re- 
moved from the group. In cases where the members of two 
groups overlapped, the two were merged into a single group, 
the linking length recalculated, and the membership adjusted 
accordingly. 

Ultimately this iterative FoF process found 1551 groups in 
the CSI catalog containing a total of 4140 confirmed galaxies 
above the mass limit, or an estimated 13000 galaxies after cor- 
recting for incompleteness (as described in the next section). 
This correction is somewhat larger than that implied by the 
mean completeness of ~ 30-40% dKelson et al.l 12012b . be- 
cause groups and clusters are relatively dense environments 
that are difficult to sample with multislit spectroscopy. Fig- 
ure Q] shows a VzJ color image of one representative group 
at z = 0.935 along with spectral energy distributions (SEDs) 
of its four confirmed members. These galaxies are circled in 
white in the image, and with this particular color scale appear 
to have consistent colors. Three of the member galaxy SEDs 
look fairly similar, with primarily evolved stellar populations, 
while the fourth is blue and star forming. This illustrates the 
high fraction of massive evolved galaxies even in relatively 
poor groups at high redshifts. 

3.2. Incompleteness and projection effects 

Even with the redshift accuracy of CSI, two effects will 
introduce serious biases if not taken into account: spurious 
groups introduced by projected alignments of galaxies, and 
group members not included in the sample due to the min- 
imum slit spacing of IMACS and yet-incomplete coverage. 
We correct for the latter effect by creating mock galaxy cat- 
alogs based on the target source density and flux distribution 
(i.e., that of all objects above the IRAC flux limit), and esti- 
mating the fraction of sources on which slits could be placed 
f(p) as a function of projected density. For a given region 
of the image, the underlying source number density is then 
approximated by n true = n obs /f(p). 

To estimate the contribution of spurious groups due to 
chance projections, we assigned random redshifts to the 
galaxies in the sample (keeping the same distribution of 
source densities) and re-ran the group finding algorithm. Not 
surprisingly, the likelihood of chance projections decreases 
with increasing group mass since a higher concentration of 
galaxies in a given region is more likely to be a "true" group 
or cluster; this contamination fraction varies from 60% at the 
low end (logM*/M Q ~ 11.5) to 10% for groups an order of 
magnitude more massive; over the full group sample the mean 
contaminant fraction is about 40%. We therefore only con- 
sider group stellar masses above log(M <r /M Q ) > 11.7, where 
the contribution from spurious groups is below 50%. 

3.3. Calculating the cumulative group stellar mass functions 

The CSI group catalog was conservatively selected from 
a galaxy sample which is complete to log(M»/M Q ) > 10.5 
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at z < 1.2 for both red and blue galaxies. Calculation of 
the GrSMFs is therefore straightforward, simply a matter of 
summing the number of groups detected in mass and red- 
shift bins and dividing by the effective survey volume (tak- 
ing into account area not covered due to gaps between masks, 
bright stars, etc.). Above the adopted mass limit of groups 
above log(M + /M Q ) = 11.7, there are 686 groups. We subdi- 
vided the sample into redshift bins z = 0.35-0.55, 0.55-0.9, 
and 0.9- 1.2. Our cumulative GrSMFs are binned such that 
the first (highest-mass) bin contains three galaxies, and the 
number of groups in each successively lower-mass bin is in- 
cremented by either one or a multiplicative factor of 1.1, 
whichever is greater; this provides a roughly constant binning 
in \ogN at low masses. Finally, to correct for projection ef- 
fects the GrSMF of the previously-described "randomized" 
(i.e. spurious) groups was calculated the same way and sub- 
tracted. 

To check the concordance of these GrSMFs with those 
found in other studies at lower redshifts and higher masses, we 
adopted two complementary samples from the lit erature: the 
Sloan Digital Sky Survey (SDSS) group catalog bv lYang et al.l 
d2007h . and massive X-ray selected clusters at z ~ and 
z > 0.35 from the ROSAT All-Sky and 400 Degree (400d) 
surveys dBurenin et al.l [20071: lyTkhlinin et all I2009al) . 400d 
includes 36 clusters at z > 0.35 and 49 at z < 0.2 detected 
in 400 deg 2 of ROSAT PSPC observations and followed up 
with Chandra to obtain total virial masses. To create a z ~ 
sample comparable to CS I, we cho s e a h ighly-restricted sub- 
set of groups from the Yang et al. (2007) "Sample II" cata- 
log with at least two galaxies of log(Mi/M Q ) > 10.5 between 
0.025 < z < 0.075 (the upper bound being the redshift where 
SDSS is complete at these masses), and calculated the cu- 
mulative group mass function. In contrast, the 400d survey 
provides an X-ray selected cluster sample given in virial, not 
stellar, masses. We therefore recomputed their mass func- 
tions in terms of Msq q using the volumes in Figure 11 of 
VikMi nhi et alj (2009a) a nd transformed the ir virial masses to 
stellar masses with the lGiodini et al.l (|2009) relations. 

If the group masses have substantial uncertainties, the ob- 
served mass function will be biased toward larger masses (due 
to the steepness of the GrSMF). This effect is small for SDSS 
and 400d, but significant for CSI which has relatively large 
mass errors. To correct for this, we convolved the SDSS mass 
function with CSI's expected stellar mass error distribution, 
and thereby estimated the shift in the CSI GrSMF as a func- 
tion of group stellar mass. All CSI GrSMFs were corrected 
accordingly, with a mass shift of ~ 0. 5 dex. 

Figure |2] shows the CSI, SDSS, and lVikhlinin et al.l d2009al) 
mass functions overplotted, all in terms of logM*/M Q . As an 
additional check, we include the H IFLUGCS X-ray s electe d 
cluster mass function from Reiprich & Bohringer (2002), 
tra nsforming the i r M20 values to M+ with relations given 
by lLagana et al.l d201 11) (which were specifically computed 
from HIFLUGCS). Poisson errors are shown on the CSI 
data points; cosmic var iance is not included, but using the 
iTrenti & Stiavellil (|2008) Cosmic Variance Calculator we es- 
timate it will add an additional 10-20% (greatest in the 0.35 < 
z < 0.55 bin and at the highest masses) to the number den- 
sity uncertainty; however, for the most part the Poisson un- 
certainties dominate. Despite the major selection differences, 
the SMFs of the z ~ optical and X-ray samples are in good 
agreement, suggesting that the adopted transformations from 
virial to stellar masses are reasonably robust. 



4. DISCUSSION 
4.1. Concordance with massive X-ray selected clusters 

The effective upper mass limit of our GrSMFs where we 
"run out" of volume are, in the realm of galaxy clusters, not 
particularly massive: at M± ~ 10 125 Mq (above which we 
only have 12 groups in our sample), in the literature they 
would typically be called "rich groups" or "poor clusters." 
The lack of rich CSI clusters is simply a result of the limited 
area covered by CSI: while 5.3 deg 2 (and the ultimate goal of 
15 deg 2 ) is very large for a z ~ 1 galaxy survey, the richest, 
most massive clusters are rare enough that much larger areas 
are needed to find significant numbers of them at low to in- 
termediate redshifts. As noted in the previous section, such 
samples are provided by large X-ray surveys covering hun- 
dred s of square degrees; h ere we have adopted the 400d sur- 
vey (VikhlininetaLL2009a) as a principal comparison sample. 

Figure [2] shows th at the cluster mass functions of 
Vikhlinin et al. (2009a) (transformed to stellar masses by the 
Giodi niet alj|2009l conversion factors) not only pick up more 
or less where CSI leaves off, but also represent a smooth con- 
tinuation of the CSI GrSMFs to a factor ~ 2 higher mass. This 
is the first demonstration of the connection between group 
and cluster SMFs at these high redshifts. While not a par- 
ticularly surprising result, it provides further evidence that 
our observing strategy and group selection methods are ro- 
bust, and that the groups found by CSI bridge a key gap be- 
tween individual galaxies and massive clusters in the distri- 
bution of dark matter halos. Although deep X-ray and IR 
observations pr ovide compelling evidence of massive clus- 
ters at z > 1 dMcCarthv et alj I2OO7I; iPapovich et ail 120101; 
iRosati et"aT]|2009t) . no homoeeneouslv-selected sample cur- 
rently exists at these redshifts; indeed, the rapid decline in 
cluster abundances at high stellar masses seen in Figure [2] (if 
it is similarly steep at z = 0.9- 1.2) suggests that such objects 
are extremely rare. The most common progenitors of low-z 
clusters must therefore lie at group masses at z > 1 . 

4.2. Hierarchical growth over the past 8 Gyr 

The combined CSI and 400d mass functions show strong 
evolution over the redshift range z = 0.35 - 1 .2, with much of 
the evolution occurring between the z ~ 1 and z ~ 0.6 redshift 
bins. Similarly, there's another significant increase in number 
density between the lowest-redshift bin of CSI and the z = 
groups and clusters. In the two lowest-redshift CSI bins there 
appears to be only marginal evolution between the mass func- 
tions; this may be a result of cosmic variance or residual mass 
uncertainties. Overall, however, the observed evolution ap- 
pears to qualitatively reflect hierarchical structure formation: 
the ranks of massive groups and clusters (logM + /M Q > 12.2) 
grow strongly over this redshift interval, while the number 
density of lower-mass groups is more constant (presumably 
because new groups are formed from below our mass limit 
as others merge into the more massive clusters). Although 
these hierarchical trends can be seen in the 400d cluster data 
alone, CSI and SDSS demonstrate that this growth continues 
to masses ~ 0.5 dex lower. 

Due to the small number of high-mass clusters, the mass 
functions in Figure [2] are shown as cumulative number densi- 
ties. To better illustrate the observed hierarchical growth, an- 
other projection of the GrSMF evolution is shown in Figure[3] 
Here the number densities of galaxies in three mass bins (two 
from CSI/SDSS and one from 400d) are shown as a func- 
tion of redshift. The difference between low- and high-mass 
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groups is striking: at log(M*/ 'Mq) =12- 12.4 the abundance 
of groups is relatively fiat, but at masses ~ 0.4 dex higher the 
number density declines far more rapidly with increasing red- 
shift. More quantiatively, we fit a power-law n ~ (1 +z) a and 
find a = -1.6 ± 0.2 for groups with 12.0 < logM < 12.4 and 
a = -4.2±0.7 at 12.4 < logM < 12.8; the difference in slopes 
is therefore significant at the 3.5a level. X-ray clusters may 
exhibit marginally faster growth (a = -5.6 ± 1.4), but due to 
the large uncertainties and potential systematics in their esti- 
mated stellar masses we cannot determine whether this steep- 
ening trend continues to logM* > 12.8. 

5. SUMMARY AND FUTURE WORK 

By employing low-resolution prism spectroscopy to obtain 
accurate redshifts for a mass-complete galaxy sample over 
5.3 deg 2 , we have assembled the most comprehensive catalog 
of 686 galaxy groups down to a stellar mass limit of 10 Mq 
and up to redshifts z ~ 1.2. The CSI GrSMFs follow a smooth 
continuation of the X-r ay selected cluster mass functions of 
Vikhl irdn et alj d2009al) when these are transformed to stel- 
lar masses. Most notably, the number density of high-mass 
groups and clusters increases more rapidly than low-mass 
groups; since the galaxies in these objects exhibit little ongo- 



ing star formation, hierarchical buildup through mergers and 
accretion must dominate their growth. 

This result illustrates the power of large prism surveys like 
CSI to effectively bridge the gap between galaxy- and cluster- 
mass halos, allowing studies of an ubiquitous galaxy envi- 
ronment. As noted before, this result is based on an early 
sample comprising about one-third of the total CSI area and 
25% of the expected spectra. While these early data pro- 
vide a qualitative picture of hierarchical growth in group stel- 
lar mass, forthcoming CSI data which attain improved com- 
pleteness over the full 15 deg 2 area will allow us to robustly 
track the stellar-mass growth of groups and clusters, and for 
the first time directly test theoretical models in the gap be- 
tween galaxy and cluster-mass dark matter halos. Addition- 
ally, deep X-ray observations (and/or stacks of existing shal- 
lower data) in the CSI fields will allow a direct comparison of 
stellar and virial masses for X-ray detected groups, providing 
an unprecedented window into the underlying halo growth. 

We thank Ryan Quadri for helpful discussions and com- 
ments on the manuscript. In addition, R.J.W. gratefully ac- 
knowledges support from the Carnegie Institution for Science 
and NSF grant AST-0707417. 
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FIG. 1 . — Left panel: Example of a group at z = 0.935 with confirmed members circled in white, shown in a 3' X 3' (1.4 Mpc on a side) region. Background 
and foreground galaxies in our spectroscopic sample are circled in brown; the large green circle is tt/2 times the group's estimated mean harmonic radius. Right 
panels: Spectral energy distributions of the four confirmed group members. The cyan histogram shows the prism data, red boxes denote broadband photometry 
and uncertainties, and the black line is the best-fit model for each galaxy. All of these objects show strong 4000A breaks, indicating their stellar populations are 
largely evolved, though one exhibits current star formation. 




l°g(M., group /M ) 

FIG. 2. — Evolution of the cumulative GrSMF from z = 0— 1. Data from CSI at low-, medium-, and high-redshif ts are show n as red, blue, and green circles 
(with Poisson error bars) respectively. The black histogram shows the z ~ mass function de rived from the Yang et al. 12007) SDSS group catalog, wh ile the 
black circles are X-ray selected z ~ clusters from HIFLUGCS ( Reiprich & Bohrineer 2002). Crosses show mass functions from the Vikhlinin et al. 1 2009a) 
400d survey converted from virial to stellar masses, with colors corresponding to the same redshift bins as the CSI and SDSS samples. All z ~ estimates of 
the mass function match up well despite the different selection techniques, and strong evolution in both group and cluster number densities is evident across the 
plotted redshift range. 
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FIG. 3. — Number density of groups in three stellar mass bins. For the two lowe r mass bins, data at z ~ and z > 0.25 are taken from SDSS and CSI 
respectively, while the most massive clusters are from the low- and high-z samples of Vikhlinin et al. (2009a), converted to stellar masses. Lines are power-law 
fits of the form n ~ (1 +z) a , with a = —1.6 ±0.2, —4.2 ±0.7, and — 5.6± 1.4 in order of increasing mass; due to potential systematic offsets between the optically- 
and X-ray-selected samples, the latter is less reliable and we show it as a dotted line. The hierarchical buildup of groups, where those with higher masses grow 
more rapidly, is evident in this figure. 



